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Abstract 

Motivated by a recent proposal on the possibility of observing a monopole 
in the band structure, and by an increasing interest on the role of Berry 
phase in spintronics, we studied the adiabatic motion of a wave packet 
of Bloch functions, under a perturbation varying slowly and incommen- 
surately to the lattice structure. We show, using only the fundamental 
principles of quantum mechanics, that the effective wave-packet dynam- 
ics is conveniently described by a set of equations of motion (EOM) for a 
semiclassical particle coupled to a nonabeltan gauge field associated with 
a geometric Berry phase. 

Our EOM can be viewed as a generalization of the standard Ehrenfest's 
theorem, and their derivation was asymptotically exact in the framework 
of linear response theory. Our analysis is entirely based on the concept 
of local Bloch bands, a good starting point for describing the adiabatic 
motion of a wave packet. One of the advantages of our approach is that 
the various types of gauge fields were classified into two categories by 
their different physical origin: (i) projection onto specific bands, (ii) time- 
dependent local Bloch basis. Using those gauge fields, we write our EOM 
in a covariant form, whereas the gauge-invariant field strength stems from 
the noncommutativity of covariant derivatives along different axes of the 
reciprocal parameter space. On the other hand, the degeneracy of Bloch 
bands makes the gauge fields nonabelian. 

For the purpose of applying our wave-packet dynamics to the analyses 
on transport phenomena in the context of Berry phase engineering, we fo- 
cused on the Hall-type and polarization currents. Our formulation turned 
out to be useful for investigating and classifying various types of topo- 
logical current on the same footing. We highlighted their symmetries, in 
particular, their behavior under time reversal (T) and space inversion (/). 
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The result of these analyses was summarized as a set of cancellation rules, 
schematically shown in Tables 1 and 2. We also introduced the concept 
of parity polarization current, which may embody the physics of orbital 
current. Together with charge/spin Hall/polarization currents, this type 
of orbital current is expected to be a potential probe for detecting and 
controling Berry phase. 

1 Introduction 

The search for a quantized magnetic monopole has a long history. ^ |2] 
Recently a group of condensed-matter physicists 0] embodied the idea of 
detecting a monopole in the band structure. |S] In crystal momentum space, 
monopoles appear as a source or a sink of the reciprocal magnetic field [HI |7] 
associated with the geometric phase of Bloch electrons. The geometric phase of 
a Bloch electron, i.e., its Berry phase, has also attracted much attention on the 
technological side, in particular, in the context of spintronics. A spin Hall effect 
has been of much theoretical concern, |S1 |H1 El HH IT^ since it may provide a 
possible efficient way to induce spin current in a semiconductor sample on which 
spintronic devices ^31 will be constructed. 

The subject studied in this paper stands at the interface between the fore- 
front of the search for a monopole and the latest technology of spintronics. We 
study the wave-packet dynamics of a Bloch electron under pertubations slowly 
varying in space and in time. We derive and analyze a set of equations of mo- 
tion (EOM) which describes the center-of-masss motion of such a wave packet 
together with its internal motion associated with its (pseudo) spin. A recip- 
rocal gauge field of geometric origin (Berry connection) appears naturally in 
such EOM. [7j Then we combine our formalism with the Boltzmann transport 
theory to describe such phenomena as spin and orbital transport. Its relevance 
to quantum charge/spin pumping |14l I15| will be also briefly discussed. 

Before going further into the detailed descritption of our project, let us 
briefly remind you what the Berry phase is, and how it has become to be widely 
recognized in the community. In his landmark paper '16?, Berry introduced it 
as a quantal phase acquired by a wave function whose Hamiltonian is subject 
to an adiabatic perturbation. The Berry connection, i.e., a gauge field appears 
as a phase of the overlap of two wave functions infinitesimally separated in the 
adiabatic parameter space. Before being formulated in such a systematic man- 
ner, the Berry phase, however, had already been recognized and discussed, for 
somewhat restricted cases though, in several independent contexts. The molecu- 
lar Aharonov-Bohm effect discussed in Ref. ^7] is nothing but a manifestation 
of Berry phase. Its relevance to band structure had also been recognized in 
limited situations, such as anomalous Hall effect (AHE), |18l I2()j as well as in 
the study of quantized Hall conductance. j2J The role of Berry phase in fer- 
roelectrics has also been of much theoretical interest. |221 Recently, the Berry 
phase in AHE has attracted a renewed attention, revealing its rich topological 
structures. |3Elil[3|23lllESll2ni The Berry phase has also been generalized 
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to a non-Abelian case. (221 

The equations of motion (EOM) for a wave packet of Bloch functions ^ is 
instrumental in all the analyses done in the paper. In order to illustrate our 
program, we begin with some details of the description of such EOM. A wave 
packet of Bloch functions is localized in the phase space around (k, x) (where 
k is a crystal momentum characterizing the Bloch function). The wave packet 
is also composed of a specific Bloch band n, whose energy dispersion relation is 
given by el^"^ (k) . The center of mass coordinates (k, x) obey a set of classical 
EOM, as the Ehrcnfcst's theorem says. In the presence of electromagnetic field 
(E, B), its motion is subject to an electric and Lorentz forces, 

-dt = ^H^w+d^^^w ' 



dx _ de^n\^ 

It ~ d\L 



(2) 



These EOM, together with the Boltzmann transport theory, describe the electro- 
magnetic response of the system. To see this point, let us express the charge 
current in terms of the momentum distribution fimction /(k) as 

where D is the dimension of coordinate space. The net current vanishes in the 
thermal equilibrium. A finite net current appears when either 

1. /(k) is deviated from its equilibrium value, or 

2. dx/dt acquires an anomalous term, i.e., an anomalous velocity. 

Case 1 corresponds obviously to the usual ohmic transport, in which the current 
is induced by a small deformation of a Fermi sphere from its thermally equili- 
brated distribution. In this case the current is, therefore, carried only by the 
electrons in the vicinity of the Fermi surface. 

The Berry phase contribution to Eq. Q corresponds to Case 2, and involves, 
in contrast to Case 1, all the electrons below the Fermi surface. This type 
of geometric current might be also dissipationless. 1281 1^ When Berry 
connection is taken into account, the classical EOM, in particular, Eq. (O is 
subject to a modification. In terms of a reciprocal magnetic field B, the EOM 
for X now reads, |7| 

where eeff(k,x, i) is an effective energy, which will be defined in more precise 
terms in Eq. (|38|l . The nature of reciprocal magnetic field will be clarified in 

Bloch function is an eigenstate of a periodic Hamiltonian such as Eq. llUi . whose 
energy spectrum forms a band structure e^\k.) defined as in Eq. Illi . 
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Sec. 2. One can observe in Eq. Q that S(k) acts quite similarly to the Lorentz 
force in the real space. S(k) encodes information on the topological nature of 
band structure, in particular, that of band crossings. 7^ Indeed, a degeneracy 
point corresponds to a momopole of i3(k), [5 which has played a crucial role in 
the understanding of anomalous Hall effect (AHE) ,61 1^^111^ . 

In this paper we study the wave-packet dynamics of Bloch electrons subject 
to a perturbation /3(x, t) varying slowly in space and time. Even though our 
treatment of /3(x, t) is in completely general terms, we can give some concrete 
examples of /3(x, t) as in Ref. 0, 

i/(p,x;/3(x,i)) =ifo(p + /3i(x,i),x + /32(x,i)) +/33(x,t). (5) 

iJ(p, x; /3 = 0) is an unperturbed Hamiltonian. The first two categories, /3i(x, t) 
and /32(x, t), are in a vectorial form, whereas /33(x, t) is a scalar. In the case of 
electro-magnetic perturbations, /92(x, i) = 0. In the following, we restrict our 
analyses, for the sake of simplicity, to the case of 

/32(x,<)=0. (6) 

A finite /32(x, i) could be relevant, e.g., for the study of deformational pertur- 
bations in a crystal. [7] 

Following the quantum mechanical motion of a wave packet localized around 
(k, x), we study its EOM focusing on the topological nature of band structure, 
and interpret them in terms of the reciprocal vector potential Aq defined in the 
(2D -I- l)-dimensional parameter space {q} — (k, x, t). This set of parameters 
{q} plays in our case the role of adiabatic parameters in the original formu- 
lation of Berry phase. JHl Our approach is entirely based on the fundamen- 
tal relations of Schrodinger qauntum mechanics, and makes no referenece to 
(i) Time-dependent variational principle (Tj, or (ii) Path- integral method using 
Wannier basis PO]. Although our approach is conceptually much simpler than 
those mentioned above, this type of analysis can be found, to our knowledge, 
only in the classical literature. |19l I^Dj We have in mind a linear response theory 
with the help of Boltzmann equation. We, therefore, restricted our analysis to 
the first order of external perturbation /3(x, t). We emphasize here that all our 
analyses are asymptotically exact in the framework of linear response theory. 

This paper is organised as follows : In Sec. 2, we first discuss the nature 
of non-Abelian gauge field, appearing in our EOM, which will be derived later 
in Sec. 4. In Sec. 3, we state and formulate unambiguously our problem, as 
well as listing all the assumptions we will make. The EOM is derived in Sec. 
4, whose possible application to Berry phase engineering is discussed in Sec. 5, 
before coming to the conclusions in Sec. 6. Some technical details are left for 
Appendices. 

2 Origin of the gauge field 

The nature of a reciprocal magnetic field B{k.) appeared in Eq. ^ lies, as 
will be further discussed in Sec. 4, in the noncommutativity of the center of 
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mass coordinates, (k, x). jJU] In more mathematical terms, S(k) is a curvature 
associated with a geometric Berry connection, i.e., a gauge field. The relation 
between such noncommutative coordinates as seen in Eqs. ()25I28|I and the MM 
in momentum space has been of much theoretical interest. From a 

more general point of view, physics in noncommutative space-time coordinates 
has been of great theoretical interest, rather in high-energy physics community, 
in particular, in the context of string and M theories. |22l In the following 
we consider, instead, the physical origins from which our gauge fields stem, and 
the mechanism how they are generated, focusing on the case of Bloch electrons 
under slowly varing perturbation /3(x, i). 

2.1 Non-Abelian gauge field, or Berry phase, encoding 
information on the band structure 

Let us consider the motion of a wave packet composed of a limited num- 
ber, say, N of degenerate bands over the whole Brillouin zone. When neither 
the time reversal symmetry nor the spatial inversion symmetry is broken, there 
always appears a two-fold degeneracy at every k-point (Kramers doublet). If 
there is no further degeneracy in the system, then = 2 in our language. In the 
following chapters, we will derive, using only the most fundamental relations of 
Schrodinger quantum mechanics, effective equations of motion (EOM) for this 
wave packet. These EOM are most conveniently interpreted in terms of non- 
Abelian gauge fields in the reciprocal space. When we derived these effective 
equations of motion, we restricted our available Hilbert space to these degener- 
ated bands. In the course of this procedure of projection onto the N bands, all 
the relevant information, about the bands integrated away, was encoded in the 
form of a gauge field, and appears in the EOM for the wave packet as a Berry 
phase. 

In order to illustrate this point, let us investigate how those gauge fields 
are expressed explicitly in terms of Bloch functions. We will see in later sec- 
tions that the concept of Bloch bands is susceptible of perturbations varying 
incommensurately to the lattice structure. As a result Bloch electrons become 
subject to a (non-Abelian) gauge field in a (2Z3+l)-dimensional parameter space 
(k, X, t), which we will call below the reciprocal space, from the view point that 
it is a generalization of the space spanned by k, the mean crystal momentum 
of the wave packet. The reciprocal vector potential takes the form of a x iV 
matrix, whose elements are given by 



where q should be understood as a general coordinate q = k^,Xv,t and fi,iy = 
1,---,Z). k,x are center of mass coordinates defined in more precise terms, 
respectively, in Eqs. ((T7|l and lfTB|) . |u„(k, x, i)) = exp(— ik • x)|0„(k, x, t)) is 
the periodic part of a local ^ Bloch state; (x + a|u„(k, x)) — (x|w„(k, x)). Inner 

^The concept of local Bloch function will be briefly introduced in Sec. 2.2 before being 
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products involving the periodic part u„(k, x, t)), mean an integration over the 
unit-cell, with a normalization (M„(k, x, t)|u„(k, x, t)) = 1. In the abelian case 
= 1, this vector potential is indeed related to the the reciprocal magnetic 
field B(k) introduced in Eq. |@J as 

In the non abelian case, the gauge invariant reciprocal field strength should be 
defined as 



9192 



— 3 A - 3 A 



92- 



where qi , q2 

last term of Eq. (jSJ can be rewriten as 



91 ^ * [-^91 ' -^92] I 

^, t. Using a trivial relation {^§^\un) + 



(8) 



|^)=0, the 
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whereas, 
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dUn 



dUn 

dqi 



(9) 



\ui) {ui\ 
would 



Comparing those two equations, one can immediately see that if X);=i 
were 1, i.e., if { ; Z = 1, • • • , N} spanned a complete basis, then J^q^^q^ 
vanish identically. This indicates the fact that the nature of our gauge field lies 
indeed in the projection of an available Hilbert space onto the relevant N bands. 
If \ui) spanned a complete basis, and no band were projected away, there would 
be no information which should be encoded in the gauge fields. Note also that 
Eq. takes the familiar form of the Berry curvature in the study of magnetic 
Bloch bands, 



2.2 Gauge field of two diff'erent origins 

The gauge field intorduccd in Eq. ((TJ has two different physical origins : 

1. Projection onto a subspace spanned by N Bloch bands, 

2. Bloch basis moving in time. 

The first point has been already discussed in Sec. 2.1, whereas the second point 
may need some explanation. In the following sections, we will study the wave- 
packet dynamics in the phase space in the presence of space and time dependent 
external pertubation, which varies incommensurately to the lattice structure. In 
order to define a crystal momentum in such a situation, we replace the spatial 
coordinate x in the perturbation /3(x, t), introduced as in Eq. |SJ|, by the center- 
of-mass coordinate x of a wave packet under consideration. This recovers the 

formulated in more precise terms in Sec. 3. 
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original lattice periodicity of the Hamiltonian, leading us to the concept of lo- 
cal Hamiltonian (Eq. (|12|l ) and its local Bloch eigenstates (Eq. (|13|) '). The 
above procedure is justified, whenever the external perturbation varies suffi- 
ciently smoothly compared with the width of the wave packet. We then expand 
the wave packet in terms of the local Bloch eigenstates, |0„(k, x(i), t)), which 
evolve as a function of time, both explicitly (through t) and implicitly (through 
x(i)). This is why our local Bloch function, or rather its periodic part, which 
has appeared in Eq. Q, depended not only on k but also x and t. Because 
of the nature of our local Bloch basis, such nontrivial gauge field structure as 
was introduced in the previous section emerges. To be precise, we had better 
distinguish between two different types of gauge field (strength) appearing in 
Eqs. 17181) : (i) .Fe z , (H) ^ and^r Although the reciprocal field strength 
introduced in Eq. ^ has various components, i.e., not only (a) Tj.^/.^, -^k^^x^ 
and J^k^t ((a)=(i)+(ii)), but also (b) Tx^,x^ and Tx^t^ the latter components (b) 
do not appear in our EOM for the wave packet, showing a clear contrast with 
Ref. 0. However, we will be working in the framework of a linear response 
theory, and within that framework our EOM turn out to be consistent, when 
TV = 1, with those of Ref. j7|. This point will be further clarified in Sec. 4.4 by 
performing a simple power counting analysis. 

We will see in detail in Sec. 4 that the two types of gauge field 

2. Ti. ^ and Ti. 

have actually slightly different origins, as well as their different physical conse- 
quences which we will discuss in Sec 5. The former, T^.^^.^, is indeed related to 
the projection of available Hilbert space onto the relevant degenerate N bands. 
It yields a finite anomalous velocity, and plays a central role in the understanding 
of AHE. It appears in the presence of magnetic Bloch bands and ferromagnetic 
backgrounds. [SI ESI 1^ • On the other hand, the latter, Ty.^^^ and T^^^ 
appear only in the presence of the time- dependent Bloch basis mentioned above. 

3 Statement of the problem 

Before discussing the EOM in the following section, let us define and formulate 
our problem here as well as listing all the assumptions we will make. We stress 
here that all the approximation which we will make are stated here, and that 
the derivation of the EOM in the following section is indeed exact under the 
assumptions made in this section. 

Let us consider the motion of a wave packet of Bloch functions under per- 
turbations slowly varying in space and time. This perturbation can be, e.g., 
external electro-magnetic field, as was the case in the study of magnetic Bloch 
bands. |21[ I5| The external perturbation /3(x, t), varing incommensurately to 
the crystal structure, breaks the translational symmetry of the imperturbed 
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Hamiltonian, 



i^o(p,x) = £- + C/(x), [/(x + a) = ;7(x). (10) 

Eigenstates of the above Hamiltonian H1U|) . i.e., Bloch bands (specified by band 
indices n) are characterized by crystal momenta k, 

i7o|0i°nk))=e("Hk)|0(")(k)). (11) 

Once the perturbation /3(x, t) is switched on, this crystal momentum k is no 
longer a good quantum number of the system. However, the typical wave length 
of the external perturbation is longer by several order of magnitudes than the 
lattice constants, in a physically relevant parameter regime of our interest. In 
that case, intermediate length scales do exist, to which our wave packet will 
belong, in which the external perturbation /3(x, t) can be regarded spatially 
constant at the zero-th order of appoximation. We are thus entitled to consider 
a wave packet, well localized in this length scale of external perturbation, which 
has also a peak sharp enough in the space of crystal momentum, moving under 
pertubations slowly varying in space and time. 

Let us now consider a wave packet, |\E'(i)), localized in the phase space, 
spanned by the real space coordinate x and the crystal momentum k, in the 
vicinity of (k, x) . For simplicity, and without losing generality, we can assume 
that the wave packet has a symmetric and smooth shape such that it has a 
well-distinguished peak at (k(i),x(t)) in the phase space, where k(i) and x(i) 
should coincide with the expectation value of k and x at a given time t. 

Our present goal is to study, as accurately as possible, the quantum mechan- 
ical motion of this wave packet, and derive the effective equations of motion for 
x and k. As will soon become clearer, an interpretation in terms of reciprocal 
gauge field (strength) uncover the nature of various physical phenomena, such as 
anomalous Hall effect (AHE), PEHIEIEHI, spin Hall effect, [EI and quantum 
charge/spin pumping. jl4l I15| 

3.1 Assumption of slowly varying perturbation /5(x, t) - 
concept of the local Hamiltonian and its local Bloch 
bands 

We consider from now on a perturbation /3(x, t) introduced in Eq. ((Sj. As far 
as the intermediate length scales discussed at the beginning of this section are 
concerned, /3(x, t) can be regarded, over the spread of our wave packet, almost 
spatially constant. We, therefore, choose, as the starting point of our analysis, a 
Hamiltonian, dubbed in Ref. |7] as a local Hamiltonian, in which x-dependence 
of /3(x, t) is replaced by x, a constant at a given time : 

ffioc =i?(p,x;/3(x,t)). (12) 

This Hioc has a very remarkable property ; at a given time t it has the same 
translational symmetry as the non-perturbed Hamiltonian Ha = H{p, x; /3 = 0), 
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i.e., in other words, i?ioc can be diagonalized by a set of local Bloch eigenstates 
|^„(k, X, t)) forming a local band e„(k, x, i), which now depends on x(t) and t : 

i/ioc \MK X, t)) = e„(k, X, t) |0„(k, X, t)) . (13) 

We are actually considering a degenerate case where e„(k, x, (n = 1, • • • , TV) 
takes the same value, which we define to be eioc(k, x,t), i.e., 

eioc(k,x,t) = ei(k,x,i) = • • • = ew(k,x,i) (14) 

We will see below that the concept of the local Hamiltonian and its assoicated 
conduction bands plays a central role in the derivation of EOM. 

3.2 Construction of a wave packet 

Superposing local Bloch functions introduced above, we now construct our 
wave packet. In the spirit of Boltzmann's transport theory, an exchange of 
energy between the electron and the environment occurs only through scattering 
events. In the following, we will investigate an adiabatic motion of this wave 
packet. This picture should be valid over the typical length scale of an adiabatic 
flight between two scattering events, i.e., over the mean free path of an electron. 
Let us now proceed step by step, making each logical step as clear as possible. 

1. Let us first focus on the real space, in which the electron wave packet is 
localized around x. Then we can compose a wave packet out of local 
Bloch functions associated with the local Hamiltonian at x : 

N 

|*(t)) dka„{k,t) |,/)„(k,x,t)). (15) 

n=l 

a„(k, t) should be normalized properly. The x-dependence of |5'(t)) is 
implicit on the left hand side of Eq. IjlSfl . which is actually due to the 
time dependent Bloch basis, |0„(k, x,<)). 

2. In order for the self-consistency, we require that our wave packet (|15ll does 
give, the correct expectation value of x, i.e.,x(t) = {xi(t), ■ ■ ■ ,XD{t)): 

x^{t)^{^{t)\x^\^{t)). (16) 

This guarantees that our wave packet yield, indeed, the center-of-mass 
position preassigned in Eq. H12|l . and that our program makes a self- 
consistent closed loop. 

Our wave packet H15() can be also regarded as a functional of a„(k, i), i.e., 
|^(t)) = |5'({a„(k, t)})), in which the coefficients a„(k, t) are chosen so that the 

•^As has been discussed at the beginning of this section, the expansion 1151 is justified, as 
far as the spread of wave packet in real space is sufficiently small compared with the typical 
length scale over which the external perturbation can be regarded almost constant. 
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self-consistency condition ((TCll should be satisfied. Eq. is, however, nothing 
but a weak constraint compared with a huge number of degrees of freedom 
allowed for a„ (k, t) . In order to specify with further precision the coefficients 
a„(k, i), we now turn our eyes to the k-space. As has been discussed at the 
beginning of this section, we can consider, in the length scale of our interest, 
a wave packet which is localized both in x and in k. We therefore require, 
in addition to Eq. H16|l . that our wave packet should also give the correct 
expectation value of fc^ : 

~k^{t) = {^{t)m^{t)). (17) 

The k dependence of is thus encoded in a„(k, t). ^ 

In the following section, we will derive the EOM for x{t) and k(t), following 
the quantum mechanical motion of the wave packet we have just prepared. In 
order to make the set of EOM self-contained, however, we need also to take 
care of the motion of internal pseudospin degrees of freedom spanned by N 
bands. For that purpose it will turn out to be convenient to separate a„(k, t) 
into its "phase" or pseudospin part, z*(k, — (zi(k, f), • • • , zjv(k, i)), and its 
"amplitude" part, p{k,t), by introducing 

a„ (k, t) = ^/p{k,t)zn (k, t) , 
p(k,i) = ^|a„(k,i) . 

n=l 

|z„P clearly represents the probability that the electron wave packet sits on 
the n-th band among the A^-fold degenerated bands. Thereby it corresponds to 
the internal degrees of freedom associated with the wave packet, such as spin 
and/or orbital, while p(k, t) is the momentum distribution function for the wave 
packet. Since we assumed that the wave packet is well- localized not only in its 
real space but also in its reciprocal space, we can assume without any loss of 
generality the following reduction formula, 

I dkf{k,t)p{k,t) = f{k{t),t), (18) 

for any sufficiently smooth function /(k, t). This prescription will be used fre- 
quently at the final stage of the derivation of EOM. 

3.3 First order perturbation theory with respect to (3{x, t): 
a hnear response theory 

The wave packet introduced above should obey the Schrodinger equation 
^ l*(t)> = H \^{t)) . (19) 

*So far, our treatment of x and k has not been symmetric. This is entirely due to the fact 
that our perturbation /3(x, t) does not depend on k. We can consider, in principle and without 
much difficulty, such a perturbation that depends on k, and perform symmetric treatment of 
X and k. However, in this paper, we restricted ourselves, for the clarity of the paper, to the 
former case. 
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As we have briefly seen in the Introduction, our eventual objective is to apply 
the EOM to the framework of the Boltzmannn transport theory, using formula 
such as Eq. ||2Jl, in order to describe phenomena including the anomalous Hall 
effect, spin Hall effect and quantum pumping, etc. For that purpose it is enough 
to consider a linear response of the system, keeping only the terms up to first 
order of 9/9(x, t)/dx and d(3{5c, t)/dt. 

In the case of the electromagnetic fields, the perturbation /3(x, t) is embodied 
by a vector potential A(x, i), and a scalar potential Ao(x); the full Hamiltonian 
reads Ho{x, p+eA(x, t)) — eAo{x.). Thereby, a linear response to applied electro- 
maganetic fields, E = —dAo/dx ~ dA/dt, B = V x A corresponds to the first 
order perturbation theory w.r.t. /?(x,t). 

We expand the Hamiltonian in powers of x — x as 

H = i/,oc + ^ E - *.)^ + ^(-. - ^.)} , (20) 

The first order term on the r.h.s. is written in a symmetrical way in order to 
keep the Hamiltonian to be hermitian. In the following, based on Eq. H20|l we 
develop a systematic perturbation theory w.r.t. /3(x,t). In this paper, we focus 
on the linear response of the system, keeping only the terms upto first order in 
the expansion. Our treatment is, therefore, self-concistent in the framework of 
linear response theory. 

4 Equations of motion 

In this section, we sketch the derivation of EOM, paying particular attention 
to, how the two different types of reciprocal field strength, introduced in Sec. 
2, appear in the EOM. Before going into the details of the derivation of EOM, 
let us remind you that there are two possible sources of Berry curvature in the 
reciprocal parameter space : 

1. projection of available Hilbert space onto the degenerated N Bloch bands, 

2. local Bloch basis changing gradually in the course of time. 

The time dependence of the local Bloch basis stems, not only from the explicit 
t dependence of the local Hamiltonian, H\oc, but also from our self-consistent 
treatment of the problem, where the local Hamiltonian depends on the center- 
of-mass position of the electron wave packet through the external perturbation. 

In many respects, our point of view is reminiscent of the standard Ehren- 
fest's theorem of quantum mechanics : the expectation value of an operator, 
such as X or p, obeys a classical EOM. We actually follow the same type of 
procedure as the derivation of the Ehrenfest's theorem, and in this sense our 
EOM can be regarded as a generalized Ehrenfest's theorem for Bloch electrons 
under perturbations varying slowly in space and time. We will come back to 
this point later. 
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4.1 Preliminaries 



We investigate, in this section, time evolution of the wave packet constructed 
in Eq. (|15|l . We are interested, not only in its motion in the phase space, but 
also in the motion of its internal spin/orbital degrees of freedom. In Sec. 5 
we will develop further analyses, from the viewpoint of transport phenomena, 
on the dynamics associated with such internal degrees of freedom. Having in 
mind applications of our formalism to those fields, we formulate our equations 
as generally as possible. More concretely, we consider the time evolution of an 
arbitary observable, O, or rather of its expectation value, 

a(i) = (*(t)|o|*(i)). 

Since we have adopted, for the sake of simplicity, the Schrodinger picture, as 
seen m Eq. Unj, the wave functions evolve in time, while observables are time- 
independent. We develop later more detailed analyses on the EOM focusing 
on the case where O = or fc^, but we consider a general observable O as 
far as possible in formulating our equations. This will make it easier to apply 
our formalism to further studies on the dynamics associated with the internal 
degrees of freedom of Bloch electron. 

Having those in mind, let us consider the expectation value of an arbitrary 
observable O, 

Oit)= J2 / dkdk'a:„(k',t)(e))_(k',k)a„(k,i), (21) 

m,n— 1 

where we have introduced an abbreviated notation for the matrix elements of 
an operator O evaluated in the restricted subspace spanned by N Bloch bands, 
i.e., 

(0)„„^(k',k) = (0„,(k',x,i)|O|0„(k,x,t)) . (22) 

Note that in this restricted Hilbcrt space not only fc^, a;^ or H but also ^ are 
considered to be an operator O. (C')„i„(k', k) is generally a N x N matrix for 
given (k',k), whereas for a given {m,n), it has, in general, off-diagonal matrix 
elements, and can be also regarded as a matrix in k-space. The presence of 
finite off-diagonal matrix elements of (C')m„(k', k), either in the k space or in 
the pseudospin space prevents some observables from commuting each other, 
thereby makes them non-abelian. 

Let us first consider two concrete examples : 

1 . Case oi O = Xfj_: The matrix elements of an observable 

(x^)_(k', k) = zS{k' - k)<5„,„-^ + <5(k' - k)(Av) (23) 

The first term is off-diagonal in k-space, when k is discrete, due to the 
k-derivative, but is diagonal w.r.t. the band index. In Eq. (|23|l . we kept 
both k' and k indices in order to emphasize the fact that this first term 
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is off-diagonal. In the following, we will omit quite frequently the k'- 
index, pretending that (x^)^^(k', k) is diagonal in k-space after S{k' — k) 
is intergrated away. On the contrary, the second term is diagonal in k- 
space, but the reciprocal vector potential Ak defined similarly to Eq. Q, 
as, 

9M„(k, x,t)" 



\ / 1 



(k,x,t) 



has off-diagonal matrix element between different bands. In the above 
equations we did not write down the explicit ^-dependence of x(i) in the 
brackets. 

2. Case of O — k^: This case is even simpler. The crystal momentum k is 
diagonal both in k and in pseudospin indices, 

(fc^)_(k',k) = <5(k'-k)<5™„fc^. (24) 

Let us further investigate the off-diagonal components of (a;p)^^^Jk', k). We 
focus here on its commutation relation in k-space. Since the first term on the 
r.h.s. of Eq. (|23|l is off-diagonal and the second term is not proproportional 
to an identity matrix, these two terms do not commute each other. One can 
indeed verify 

{xf,),{x^) (k) ^iTk^k,,, (25) 

J mn 

where [A, B] = AB — BA is a standard commutator of two N x N matrices 
A,B. Thus the noncommutativity of (a;^)^^^(k) turns out to be the origin of 
the emergence of Tk^k^- Another important remark on the Eq. (|23|l is that it 
leads us to introduce naturally the concept of covariant derivative in momentum 
space 20], defined as 



mn 



In terms of the covariant derivative ( Vt ) thus introduced, the matrix ele- 
ments (a;p)^^^(k) can be rewritten simply as 

(-.)mn(k)-<V.,)_. (27) 

The commutator between two covariant derivatives along different axes is di- 
rectly related a non-abelian Berry curvature ; 



mn 



'i^k.k^- (28) 



In geometric terms, Eq. H28II can be interpreted in such a way that two parallel 
transports along different axes on a curved surface generally do not commute 
each other. 
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4.2 To derive the EOM 



Let us now consider the time derivative of the expectation value, 0(t). Expand- 
ing in terms of the local Bloch functions as Eq. (|15|l . one can classify the 
time derivative of 0{t) into three parts : 



dOjt) 
dt 



E 

cr,o-" 



dt \ 



,(k,k")a."(k",0 



+ T.J dk'dk"a:,(k',<)(|(0)^,^„(k',k"))a..(k",t) 

+ J2j dk'dkal,{k',t){0)^,^{k',k)^^^. (29) 

a' (T 

We have in mind that the operator O is either x^, fc^j or some other observables. 
In the case of standard Ehrenfest's theorem, 

the second term of Eq. H29() docs not exist, since the matrix element, (O) (k', k) 
is time-dependent only when the local Bloch basis evolves in time. The first and 
the third terms, i.e., the change of expansion coefficients ag-(k,t) yields a com- 
mutator, [iJ, O]. They contain, however, also a Berry connection contribution, 
which, together with the second term, produce a new type of contribution, which 
we will call il^"* in Eq. (|33() . The first term of Eq. H33|l . is a generalizaton 
of the standard Ehrenfest's commutator [H,0], which induces, when O = a;^, 

(2) 

Tj. j.Jvci the EOM. On the other hand, the second term, fix^ , can be rewritten 



in terms of Tz. ^ and Tr. , . 



In order to rewrite the Eq. H29|l in terms of fi^' and let us first look 
into the following relation. 



g<(k,t) 
dt 



N 



= ^a:„(k,t) 



idt)r 



+ i(H) (k) 

\ / ran ^ ' 



(30) 



The first term is a Berry connection contribution. As is clear when it is written 
more precisely as 



(1) 

\ OT I mil 



dt 



■ti„(k,x(t),i) — u„(k,x(t),i) 



(31) 



it emerged as a result of the time evolution of the local Bloch basis. On the 
other hand, the second term of Eq. H30|l yields the commutation [H, 0]„m in 
Eq. H33I) . Note also that the derivative d/dt in Eq. H3U|) picks up both the 
explicit and implicit t-dependence. Correspondingly, one can also rewrite the 
first term of Eq. (jSOj) using two types of the gauge field introduced in Sec. 2, 
i.e., 

dx,. 



(I) 



dt 



At 



(32) 
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Our next objective is to calculate the time derivative of an operator such as 
k, X and express them in such a way that their interpretation in terms of the 
reciprocal field strength will become as easy as possible. For that purpose, we 
rearrange the terms in Eq. (|29|l into two parts, ril^' and fil^^ as, ^ 



(1) 

o 



n 



^0 1 



N 



m,n— 1 
N 



J2 dk^dk2a*^{ki) {H),{0) (ki,k2)a„(k2 



^ / dkidk2a'^{kx) 



(O) (ki,k2) 



+ (^(0)m„(ki,k2)) la„(k2) 



(33) 



In Eq. H33() we did not write down explicitly, for the sake of simplicity, the 
dependence on x and the t in |(/)„(k)) = |0„(k, x(i), i)). 

As has been announced in advance, i}^^ is a generalization (or, rather a 
restricted version) of the standard Ehrenfest's commutator, coming exclusively 
from the first and third terms of Eq. H29|l . whereas fi^'' is a new type of 
contribution, which is a collection of Berry curvature terms from all the three 
parts of Eq. H29|) . Not only have they different origins, but also are they 
susceptible of different physical interpretations in terms of the reciprocal field 
strength. We will see in Sec. 4.2 that in the particular case of O = x^, the two 



contributions, flx''^ and are related actually to the two different parts of 
the gauge field introduced in Sec. 2, i.e., (i) J^k^k^, and (ii) J^k^x^, ^k^t- 

Having in mind what has been stated above, we can now derive the EOM for 
x^{t) and kfi{t). Let us first consider the case of O = fc^. As seen in Eq. (|^ . 
the momentum operator is not only diagonal in k coordinates and band 
indices, but also its matrix element is time- independent. Thus only the first 
term Jl^,^'' contributes to its EOM . Futhermore, among various matrix elements 
of the Hamiltonian given in H36I37|I . only those terms which contain off-diagonal 
matrix elements w.r.t. k indices contribute to its commutator with k^. As a 
result, its EOM turns out to be simplified as. 



'^^-(')=-/dkp(k,t)-^^'°^(^'*'*) 



dt 



5eioc(k,x, t) 



(34) 



has only to substitute literally Eqs. I3UI into the 



^In order to obtain Eq. 1331 . 
expression for dO/dt in Eqs. J^, and rename the dummy variables in the following way 



k' 



ki, k 
m, a 



■ k2,k - 
n,cr - 
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In the second equality, we replaced k in the integrant by its mean value, following 
the prescription given in Eq. 10 . This is nothing but the standard EOM for 
the momentum of the electron wave packet shown in Eq.JQ). 

As for the position operator , Eq. (|23|) contains both the k-derivative and 
time-dependent matrix elements between different band indices. As a result, the 
EOM for the real space coordinate is subject to a drastic change in comparison 
with Eq. |5J. In Sec. 2, we classified the reciprocal fields into two categories, 
i.e., (i) Tk^k^, and (ii) Tk^t and J^k^x^^§f- We will see in the next section 
that the decomposition clearly demonstrates why we classified them in 
that way. We have studied on a very general basis in this section that the 
two components, and are structurally well distinguishable, and have 
completely different nature. We will see more specifically in the next section 
that r^i^^' and ili?^ are realted respectively to the reciprocal fields (i) and (ii). 
Thus different origins of two types of reciprocal fields will be uncovered. The 
fact that the classification of reciprocal fields discussed in Sec. 2 can be done 
explicitly and unambiguously as the decomposition (|33l) . is actually one of the 
main advantages of our approach. Let us now turn to a close inspection of the 
nature of fl^^ and Q^x} . 

4.3 Nature of fiW and r^l^) 

In this section let us further analyze the nature of decopmposition (|33|) . 
focusing on the case of O = cc^. x(i) is given by Eq. lfTC|l together with Eq. 

The first term of Eq. i}^^ in the present case, has particularly a 

familiar form, which often appears in the context of the Ehrenfest's theorem 

I {^\x^\^) ^ I [H, x,]\^P) = ^^1^1^) . (35) 

The similarity between and Eq. (|55|l becomes clearer, when one expands the 
wave packet lip) in terms of a complete set of bases as \4>) = \4'a)- 
The difference is that the set of bases used in the expansion was complete in Eq. 
H35|) , whereas it was restricted to N Block bands in fl^x^ . This constraint is the 
origin of non- vanishing field strengths. 

Let us now proceed to rewrite Jli^^^ in terms of the reciprocal field strength, 
J^k fc„, defined in ©. The matrix elements of the Hamiltonian, i.e., Eq. H2U|) in 
the restricted subspace, spanned by N Bloch bands are calculated to be 

(i7)_(k',k) = (0,„(k',x,f)|iJ|0„(k,x,t)) =<5(k'-k)(i/)_(k) (36) 




where ejoc = eioc(k,x,i) is a degenerate eigenvalue of the local Hamiltonian ^. 
We also introduced a renormalized energy, eeff(k,x,i), which takes the form of 

•^See Eqs. 1121141 . Recall also that 

(Hioc) (k',k) = (</.„(k',X,t)|i?locUn(k,X,t)) =5(k' -k)5™„eioc(k,X,t) 
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a iV by TV matrix whose (m, n)-components are given by 



e5^„(k, X, t) = eioc(k, X, t)dmn + AemniK X, t). 

Its off-diagonal matrix elements are due to correction terms, 



- .^ _ i I gMm(k,x, t) 

i I duyn{]&,yi,t) 
2 



(tj /, - 5u„(k,x,<) 

[H\oc - eioc(k,x,t) ' 

(i^loc — eioc(k, X, t) 



5u„(k,x,t) 



dk,, 



(38) 



where the summation over ^ = 1 , • • • , I? was assumed implicitly. Using the 
matrix elements given in Eqs. (|36I37|I one can rewrite Jli^^^ in the following way, 



fill) = 



N 

E 

m,n— 1 
D 



rfkp(k, i)z;„ (k, t) I I^Vfc^ , Ecff (k, x, t) 
5eioc(k,x, i) 



z„(k,t) + Af]«, 



(39) 



Apart from the energy correction Ae, the first term of ^l^x^ is nothing but a 
standard velocity term, i.e., the first term in the r.h.s. of Eq. A remark 

worth mentioning here is that the covariant derivative in the commutator plays 
a central role in ensuring the SU{N) gauge invariance of final results, which 
we will see later. Afii^^^ is a irrelevant term ^ which vanishes with the help of 
prescription introduced in Eq. I|18|l . For a later convenience, let us introduce 
the following abbreviated vector notation for Zn(k(<), i) : 



B(t)t= (z*(k(t),t),...,z^(k(t),t)), 



(40) 



where we always have in mind the prescription 118|l . Using this notation, one 
can further rewrite Eq. as 



JV 

E 

m.n—l 



Vt ,eeff(k) 



l/=l 



9eioc(k) 
I dxi, 



Zn{t), (41) 



where x-dependence of eioc(k) is not written explicitly. 

In contrast to the second term of Eq. (|33|) . fl^x,! in the present case, 

would not have existed, unless the local Bloch basis had evolved in time. How- 
ever, in a general situation described by a time-dependent Bloch basis, there is 

(2) 

no reason to believe that Jli,^ should vanish. Indeed, we will give you in Sec. 

is proportional to an identity in the pseudospin space. 
'^Details are given in Appendix A. 

*Its explicit form is given in Eq. 11021 in Appendix A. 
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5 some concrete examples where a finite contribution from flxj^ plays a crucial 
role in determining the physical properties of the system. Using the Eqs. (p^ . 

(2) 

one can easily verify that flxj^ are related to the second category of reciprocal 
fields, i.e., Tk^t and Tk^,x^^^- ^x'} can be rewritten as 

= -zWt|f]^,^^,^^+^,^Az(t), (42) 

where the summation over v — 1,---,D was omitted in the first line. The 
decomposition together with Eqs. (|^ ESJ gives a complete physical justi- 
fication of the classification of Tq-^q^ done in Sec. 2. In other approaches |7I \'M)\ 
the two types of reciprocal fields appear in an indistinguishable manner, and 
two different origins of reciprocal gauge field studied in this paper remain to be 
hidden. 



4.4 SU{N) gauge invariance 

We have successfully related the two contributions to ^x(t) in the decom- 
position (|^ . i.e., il^x^ and respectively, to two types of gauge invariant 
reciprocal fields, (i) J^k^,k„, and (ii) J^k^t and J^k^,x^^ ■ Together with Eq. 
this allows us to rewrite our EOM for x(i) and k(i) as 



dx n _+ 
— — = z' 



dt 

dkfj. _ 9£ioc(k, x,t) 
dt dxfj_ 



- ^k.-k.^ - -^S.s.^ - -^S^tl^, (43) 

(44) 



Repeated indices v should be summed over v — 1, - ■ ■ ,D. The effective energy 
eeff is related to the local eioc as Eq. In Eqs. 143I44|) . €cs and eioc are 

functions of k,x, i.e., CcS — ecff(k, x, i), eioc — eioc(k, x, t). Below let us verify 
explicitly that Eqs. H43() are indeed SU{N) gauge invariant. In order to obtain 
a complete set of EOM, however, we still need to know an EOM for z{t) defined 
in (|40|l . The details of its derivation is given in Appendix B, and the result is, 



where repeated indices /i should be summed over /i = l,---,Af, and again, 
Ccff — Ccff (k, X, i), eioc — eioc(k, x, t). Eqs. H43I44I45II consitutc the central result 
of this paper together with their interpretation in terms of the field strength 
•^9192) which will be further discussed in Sec. 5. 

Before turning to the discussion on gauge invariance, let us make a few 
comments on Eq. (|45|) . On the r.h.s. the diagonal part of ecfr, i.e., eioc(k)l, 
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tpgethcr with the second term simply give rise to a usual U{1) phase factor 
associated with an effective energy, eioc(k, x, i) — ^"^'"g^,*^'^'*^ . On the 

other hand, Ae(k) generally has off-diagonal matrix elements between different 
bands and thereby yields a nontrivial SU{N) phase factor, which corresponds 
to the precession of the spin and/or orbital associated with the wave packet. 
The remaining terms of Eq. I|45|) represents a Berry- Wilczek-Zee phase |77] 
originating from the adibatic motion of the wave packet. The first two terms 
are due to its motion in (k,x)-space. In the Abelian case {N — 1), the EOM for 
Tc{t) and k(t), i.e., Eqs. H43I44(I . are independent of the motion of phase degree 
of freedom, z(t), whereas z{t) acquires a quantal phase due to the evolution of 

X and k : exp [i J dtl^^^A^^^ + -^Ax^ +At)\ where the summation over /i was 
omitted. This is analogous to the Berry quantal phase [lb|. 

Since the A^-fold Bloch states are energetically degenerate over the whole 
Brillouin zone, these EOM should be independent of the choice of N Bloch 
bases and be invariant under the following gauge transformation : 

N 

|u„(k,x, i)) = ^ |u,„(k,x,t))5r„i„(k,x,t) 

m— 1 

z(k,t) = .g"i(k,x,i)z(k,i). (46) 

Here |u„(k,x, i)) and z(k,i) are transformed inversely to each other, making 
the l.h.s. of Eq. H15|l invariant. The gauge field and the field strength associated 
with it are transformed in the following way. 



-1 dg 



where = k, x, t. Concomitantly, the covariant derivative, defined generally 
for this as 

V,. = |--^A., (47) 
obeys the following transformation rule: 

V,. = g-'^q^g- (48) 

Futhermore, one can also check that the N x N matrix ecff (k) is transformed as 

ecff(k) = .g^^ecff (k)g. 

Using the above transformation rules, one can indeed verify that our EOM H43L 
144(1 are invariant under SU{N) gauge transformation H46() . 
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4.5 Abelian case : comparison with other approaches 

In the Abelian case : = 1, the above equations of motion 143I44|) reduces 



to 



dx^ _ i9ecff(k,x,i) ^ _ t ^ _ -r- (aq\ 



dt dk^ dt dt 



dk^ _ deioc(k,x,t) 
dt dx^ 



(50) 



EOM similar to Eqs. H49I50|I have been derived, to our knowledge, twice, using 
either 

1. Time-dependent variational principle |2] or, 

2. Path-integral method using Wannier basis. ^ 

If we compare Eqs. (2.19) of Ref. |7] and our Eqs. (|49I50|I . it can be observed 
that three terms, 

^s.s.^+^,,fe.^+^s,t (51) 

are lacking on the right hand side of (|50|l . However, one can easily check by a 
simple power counting that these terms appear only at orders higher than 2 in 
the perturbation series w.r.t. or x — x. Let us briefly illustrate this point. 
Since a subscript x implies a derivative w.r.t. x, which is always accompanied 
by X — X, it increases the power by one. It is also the case for the subscript 
t. Therefore, the first and the last terms of H50|) turs out immediately to be 
at least of the second order of (3. In the second line, i.e., in Eq. H50|l . since 
9eioc(k, X, t)/dXfi is also of the first order w.r.t. one can verify that dk^^/dt is 
at least of the first order w.r.t. f3, even if (|51|l is added to it. Taking everything 
into account, one can conclude that the lacking terms H51|) are at least of the 
second order w.r.t. /3 or x — x. 

Another difference between Eqs. (|49I50|I and Eqs. (2.19) of Ref. jT] is that 
in our EOM for dk/dt, the derivative d/dx^ applies to eioc(k, x,f) and not to 
eeft(k, X, t) as in Ref. 7*. In our formalism, as is clearly shown in Eq. 1)34(1 . there 
is no room for Ae(k, x, i) to enter the expression (|5U|I . Nevertheless, repeating 
the same type argument, i.e., the power counting for Ae(k, x, t), one can confirm 
that the contribution from Ae(k, x, t) is not physically relevant at the first order 
of /3 or of X — X. 

We have not only developed a systematic perturbation theory w.r.t. /3 or 
X — X, but also we make no approximation apart from the assumptions stated 
in Sec. 3. Our calculation must be, therefore, exact at the first order of per- 
turbation theory. Since possible discrepancies start only at the second order in 



^We had some difficulty to jusitify the use of Wannier basis used in Ref. \?iVi\ as a complete 
basis necessary in the path integral formalism. This is closely related to the arbitrariness of 



Wannier function discussed extensively in Ref. 1341 
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the perturbation series, our result, Eq. (|49I50|I is not inconsistent with that 
ofRef. 

5 Discussion : Berry phase engineering 

The gauge invariant EOM 143I44II have been successfuUy derived in the 
previous section. The decomposition (|33|l uncovered the orgin of two different 
types of reciprocal fields introduced in Sec. 2. In this section we discuss some 
physical consequences of Sec. 4 in the context of Berry phase engineering. 

In the Introduction, we argued that a finite net charge current could be in- 
duced by the U{1) Berry phase correction to the semiclassical EOM (|1I2|) . This 
finite charge current is actually carried by all the electrons below the Fermi sur- 
face, i.e., by the electrons in the ground state. Generalizing this U{1) argument 
to the non-Abelian case, we will discuss in this section how the various types of 
non-abelian field strength appearing in our EOM are related to concrete physical 
realizations, mainly focusing on the SU{2) case. This opens a new possibility 
of manipulating the ground state electronic wave function by controlling Berry 
phase, which is sometimes called. Berry phase engineering. 

After introducing some terminologies and fixing notations, we will focus on 
two topics. In Sec. 5.2, we will see that Tk^k^ is related to the physics of Hall 
type current. We first observe that the charge Hall current can be described by 
a trace of non-abelian field strength J^kj^k^ , while this current vanishes whenever 
the system is time-reversally (T-) invariant. The charge Hall current carried by a 
(k, t) Bloch electron and that of the (— k, |) electron precisely cancel each other. 
Based on this observation, we then propose two physical situations in which the 
difficulty of the cancellation of Hall type charge current will be overcome and the 
reciprocal magnetic field manifests itself experimentally as a Hall type current. 
The two situations arc : 

1. anomalous (charge) Hall current observed in systems with broken reversal 
symmetry, i.e., in ferromagnets, p i^Oll^ 

2. spin Hall current in time reversally symmetric systems. [51I^ I10[I11[I12| 

In Sec. 5.3, we will argue that J-k^^t is directly related to various types 
of polarization currents, currents induced in insulators under time-dependent 
perturbations. Similarly to the case of Hall type charge current, which will be 
discussed in Sec. 5.2, we first observe that in such systems that are symmetric 
under spatial inversion, the polarization electric/spin current actually vanishes 
due to a cancellation associated with the inversion symmetry of the system. 
Then, in order to overcome this difficulty we propose, in parallel with Sec. 5.2, 
two physical systems in which the problem of cancellation will be resolved, and 

^"^The calculation done in Ref. Q is also first order, in particular, their Hamiltonian, Eqs. 
(2.1), (2.14) and (2.15), is first order, but they kept all the possible Berry phase contribution 
without making further consideration of power counting, whereas we omitted systematically 
higher order terms. 
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the gauge invariant reciprocal field strength appears explicitly in a macroscopic 
physcial quantity, i.e., as a polarization current: 

1. If the inversion symmetry is broken externally or spontaneously, the abelian 
Tk^t gives rise to a relevant contribution to the polarization electric/spin 
current. 

2. Even in systems symmetric under spatial inversion, non-abelian Tk^t may 
have a chance to manifest itself as a polarization orbital current, if that 
orbital degree of freedom changes its sign under the spatial inversion. 

The analogy and correspondence between the Hall type and polarization cur- 
rents are summarized in Tables 1 and 2. 



5.1 Preliminaries 

Before further discussing Berry phase transport, we first introduce some ter- 
minologies as well as giving an unambiguous definition to spin/orbital currents. 

In order to illustrate our point, let us first consider a spin current. We 
have naturally in mind that there are different points of view |1(J[I11[[T^ on the 
definition of spin current operator, J^^. The difficulty of defininig a spin current 
stems simply from the fact that the bare spin Sa is generally not a conserved 
quantity due to spin-orbit interaction, i.e., the continuity eqation, dSa/dt + 
^^^-^ dJ^^/dXfj, = 0, is not satisfied. Thereby the Noether's theorem does not 
apply. If one focuses on the time derivative of a local spin Sa (x, t) in the general 
case of non-conserved spin, one could observe that there are two contributions 
to it of physically different nature, i.e., contributions from (i) a local spin curent, 
(ii) a local precession of spin. The former is the one of our interest, and the latter 
is related to the non-conservation of spin. Unfortunately, there is no systematic 
prescription for distinguishing between those two contributions. 

We can still define on quite general ground a current operator Jj associ- 
ated with an internal degree of freedom X as the time derivative of a spatial 
polarization of X as. 



'i = 



dt ' 

M 



(52) 



= JdJ2 li^^^'^+^'j^i}- (53) 

i=i 

L and AI denote the system size and the total number of electrons, respectively. 
The subscripts j attributed to I and x specify each electron. In the case of 
a spin current, the operator I should be replaced by a usual spin operator, 
J = Sa- The main reason why we have defined the spin/orbital current as 
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Eq. (|52I53|) is that that is the one which could be diretcly observable. For 
example, an increasing (-Ps„)^ = Pf_i{Sa) indicates that extra up (down)-spin 
electrons with Sa = +1/2 (—1/2) acculmulate in one (the other) end of a system 
with Xf^ — L/2 (— _L/2), which could be experimentally detected by some optical 
probes. 

We will also discuss orbital currents. A Bloch electron has, in addition to 
the spin degree of freedom, orbital degrees of freedom 11 in multiband systems. 
These orbital degrees of freedom describe the charge distribution in the unit 
cell, and hence, 

n(x,p) =n(x + a,p). (55) 

In Sec. 5.3 we focus on an orbital operator 11 which behaves quite contrastingly 
to the spin under time reversal and spatial inversion, i.e.. 



n(x,p) 
n(x,p) 



-n(-x,-p), 
n*(x,p). 



(56) 
(57) 



In contrast to the spin operator S, the parity operator 11 reverses its sign under 
spatial inversion, while invariant under time reversal. Accordingly, we dub this 
oribtal operator as a parity operator. Various transformation properties of 
the operators X — Sa,Tl are summarized in Table 1 . 

The spin/orbital current associated with an operator 2 was introduced in 
Eqs. (|52I53|) . We calculate in Sees. 5.2 and 5.3 those spin/parity currents 
carried by the ground state. In order to make later discussions clearer and 
physically more appealing, we make the following assumptions: 

1. We keep only their matrix elements in the restricted subspace J\f spanned 
by TV-fold degenerate bands ; If to e A/" and I ^ Af, then 



(0„,(k)|5„|</)Kk)) = 0, 

,(k) n(x,p) M^)\ 



(58) 
(59) 



Otherwise, we could have defined it also as 

M 



1 V 1 J dx cZx 

This type of definition is convenient for the application of Kubo formula. [121 1111 1291 

^•^The periodicity of 11, i.e., Eq. 1551 excludes its off-diagonal matrix elements in k-space: 

(n)„„(k',k) = 5(k' - k)(n),„„{k) 



(n>„„(k) 
_ (2^)^ 



n 



0n(k,i) 



Junit coll 



cix ^nm(k, t) x^n^ 



M„(k, t) ). 



(54) 



^^This orbital operator is different from the angular momentum operator, lj = x X p, which 
we might also call an orbital operator, oj reverses its sign under time reversal, while remains 
invariant under spatial inversion. 
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As will be seen in Sees. 5.2 and 5.3, these approximations make the follow- 
ing discussions considerably simpler, i.e., not only a charge current but also 
spin/parity currents become related simply to the non-abelian field strength, 
^k,,K and J^fe,,t. 



5.2 ^k^K induces Hall type currents: AHE and spin Hall 
effect 

The field strength Tk^k^ describes various kinds of spontaneous Hall currents 
carried by the ground state. In order to demonstrate it we expose our system 
under a uniform electric field E; we consider below the following Hamiltonian, 

D 

i7(p,x,/3(x,<)) =i/o(P,x)+e^£;^x^. (60) 

We see below that the Hall-type current can be expressed essentially as a trace 
of TTk^k^ in the TV-fold degenerate pseudospin space. The cancellation or the 
survival of such Hall-type topological current is essentially related to the trans- 
formation properties under time reversal (T) of the system. In systems with 
broken T symmetry a finite charge/mass Hall current generally appear. In two 
spatial dimension 13 = 2, this situation is often described in terms of Chern- 
Simons gauge field, accounting for the quantized Hall conductance, |35| as well 
as fractional charge and statistics. Chern-Simons terms also appear in 

electrically neutral systems. |37] Quantization of spin Hall conductance in un- 
conventional superfluids has also been studied in this context. j38| 



5.2.1 Charge Hall current 

Let us first see that the trace of the non-abelian field strength Tk^k^ describes 
a spontaneous charge Hall conductivity, i.e., an anomalous Hall conductivity, 
ini 123 IMl im In terms of Eqs. H52I53I) . we are considering the case of X = 1. 
Applying the EOM, Eqs. (|43I44|I . to the present case, we consider a Bloch 
electron, in the j-th pseudospin state and with a crystal momentum k. The 
SU (N) gauge invariant EOM for this electron under a uniform electric field E 
now reads, 

da;^-'^ 9eioc(k) ^ (i)f t /r\ (i) '^^•^ fat\ 

IF - -dk;r-2:-''''^^M^)-''' (61) 



~dt 



= -eE^. (62) 



When the crystal momentum k is located below the Fermi surface, all those 
pseudospin states are completely occupied and contribute to the charge current. 

^''The j-th eiegenstate is a linear combination of different pseudospin states m = 1, ■ ■ ■ , N , 
i.e., X/m— 1 l'^m(k)). These eigenstates are chosen to be orthogonal to each other: 
z(')tzW) = Si^ for j,j = 1,---,7V. 
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In order to calculate the charge Hall current, we need 



= -iVe^^|M - Tr (k)] i?.. (63) 

Then, by integrating these contributions over filled k points, we obtain the total 
current carried by all the electrons below the Fermi energy ep, 

- E / , ^ 7^ Tr E.. (64) 

Here we assumed for the sake of simplicity that the non-perturbed Hamiltonian 
-ffo(P:X) is also invariant under spatial inversion (/). As far as the Hall type 
current is concerned, however, this / symmetry plays only a minor role. Eq. 
H64|l takes indeed the form of a Hall type current reflecting the antisymmetry of 
J^k K'- fc„ — —Pk^k ■ The first term of Eq. (|63|) did not contribute to Eq. 
H64|) due to a cancellation associated with the / symmetry. On the contrary, 
Eq. H64|) remains finite irrespective of the /-invariance of Hq. 

Unfortunately, the charge Hall current obtained in Eq. I|64() vanishes when- 
ever the system is T-invariant. Let us see this point more explicitly. We con- 
sider an unperturbed Hamiltonian _ffo(x, p) which is invariant under T. The 
T-invariance relates its A^-fold degenerate Bloch functions at k with those at 
— k upto a certain SU{N) gauge transformations .g(k), 

N 

In the language of field strength, this reduces to, 

.Ffc,,fe„(-k)*--g-i.F,^fe„(k)5, (66) 

where the superscript t represents a transposed matrix, i.e., = {^)nm- 

One can verify this using Eqs. H7I8I65|I . Consequently, the charge current carried 
by a Bloch electron at k cancels with that of — k electron, i.e., 

D D 

J2 Tr (k)] E,^-J2 Tr (-k)] E,. (67) 

1^=1 v=l 

Eqs. (|64I67|) indicate the absence of spontaneous Hall current in T-invariant 
systems. 

5.2.2 Spin Hall current 

We are thus led to investigate the spin current, expecting that the spin 
Hall current remains finite even in T-invariant systems. |H1 El EB ^] The 
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underlying idea is simply that such a sign change as seen in Eq. may 
be compensated by that of spin operator under the operation of time reversal. 
The spin current has been defined in Eq. H52|l , in a more general context for an 
arbitrary internal degree of freedom X. Let us first observe that the T-invariance 
of i?o (x, p) is instrumental for this compensation. The total spin carried by the 
Bloch electrons at k has the same absolute value and opposite sign of that of 
the Bloch electrons at — k, 

Tr[(5„)(k)] =-Tr[(5„)(-k)], (69) 

where (5*^) (k) is a A'^ x matrix, whose (m, n)-components are given by 

(5„)_(k) = (0„(k)|5„|0„(k)) 



(27r 



)° [ 1 

— / dyi (w„i(k)|x,a)-(cra)^j^(x,6|u„(k)). (70) 

11 J unit coll „ I ^ 



Let us now consider the spin Hall current defined as (|52I53|I with I being I — Sa, 
the usual spin operator. In order to evaluate a spin Hall current, the EOM 
I|61I62|) used for the calculation of charge Hall current are no longer sufficient. 
We need instead to derive EOM for an observable, 

0% = O^iS^) = l{S^x, + x^S^). (71) 

Following the same type of procedure as the derivation of Eqs. H61I62|) . we 
perform, in particular, the decomposition H33|l . As was also the case in Eqs. 
H61I62|) . the second contribution to H33|l . i.e., fi^'' vanishes in the present case. 
^6 The EOM reads, 

d&fia _ . \ ^ / Ji,ji,/ * n, +^ / Tj\ / l''^S 



dt 

m.n— 1 



/ dkdk'a:jk,t)[(i7),(0^„)l (k,k')a„(k',t) 
/ dka;. (i7o) + e J2 E^{x^),{Ol^) a„. (72) 

J L ^— ' J mn 



N ^ D 
m,n— 1 fi—1 



The assumption (|58|) allows us to rewrite Eq. H72|) as 



— - ^ , dkp(k)z;„(k,i)[(5„)„„(k) 

m,n— 1 

D 

+ i E { ( (^") (k)) ^^E, + h.c. }] z„(k, t). (73) 

1^=1 ™" 



^^To be more specific, T-invariance relates the matrix of Sa at k and at — k through a 
SU {N) gauge transformation g, 

(5„>(-k)* = -<;-i(5c>(k)3. (68) 

Eq. I68i can be verified explicitly using Eq. I65i . 

^^As long as the electric field E in Eq. l6Ui is uniform, the local Bloch function defined in 
Eq. I12i has no dependence on x and t. As a result, Q^? vanishes. 
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The details of the derivation of Eq. (|75|l is given in Appendix C. We then 
apply the prescription 1)18(1 . replacing k in the integrant by its mean value k. 
The contribution to the spin Hall current by an electron occupying the j-th 
pseudospin state at k is, therfore, 



dt 



+ 



E 



0) 



(74) 



Finally we take the summation over N pseudospin states and over filled k points 
to find, 



' flOt 



dk 



^ d(¥^^^ 



E 



€ioc(k)<eF 



dk 



dt 



Tr 



{S^){k)TkM^) E.- (75) 



As was the case in Eq. (|64() . we also assumed in Eq. (|75|) that iJo is invariant 
under /, so that the first term of Eq. H74|l should not appear in Eq. (|75f) . i.e.. 



eioc(k)<eF 



dk 

{2nY 



■Tr 



{So) (k) 



5eioc(k) 



0. 



(76) 



Contrary to this normal part, the spin Hall current associated with the anoma- 
lous velocity, given by Eq. 1751) . has a possibility to be finite irrespective of 
the / and T symmetries. When the system is invariant under either T or /, 
the spin current carried by a Bloch electron at k and that of — k give the same 
contribution; 



Tr {Sc){k)TkM^) =Tr {So) {~\<i)TkM'~^. 



(79) 



Under T symmetry Eq. (|7^ is a consequence of Eqs. Ht)()l()8|) . Eqs. H75I79|I 
confirm that our expectation that the spin Hall current is robust against T 
symmetry was indeed the case. A set of cancellation rules for the Hall type 
currents are established in Table 2. 



^''^In the /-invariant case, instead of Eq. 1651 . Eq. 1851 holds, i.e., the 7-invariance of Hq 
relates — x,a|Mi(~k)) with (^x, r | Mj (k)) up to a SU(N) gauge degree of freedom h. This 
reduces in terms of spin and field strength to 

(S^X-k) = /i-i(S„)(k)h, (77) 
•5^ft^fc„(-k) = /i-'.Ffc^fc„(k)h. (78) 
Eq. 1771 justifies 1761 . whereas multiplying Eqs. I77i and I78i . one finds immediately Eq. 
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5.3 Tk^,t induces a polarization current: parity polariza- 
tion current and quantum spin pump 

We have seen in the previous section that Tk^k^ is related to Hall type 
currents associated with the internal degrees of freedom such as charge and 
spin, by applying a uniform electric field to the system. Here we argue that J^k^t 
describes various kinds of polarization current. More specifically, we consider a 
situation where a band insulator is subject to a time-dependent perturbation 
which does not break the periodicity of the underlying crystal, i.e., we consider 
a Hamiltonian, 

i/(x,p;/3(t)) -i/(x + a,p;/3(0). (80) 

Since the perturbation P{t) does not depend on x, the local Hamiltonian defined 
in Eq. (|12|l reduces simply to 

i?ioc(x,p,<) = iJ(x,p;/3(t)). 

Depending on the perturbation P{t), the ground state wave function of the lo- 
cal Hamiltonian, iJioc(x, p,t) also evolves temporally. Since an electronic wave 
function for the ground state describes spatial distributions of charge, spin and 
orbital, its evolution in general induces various kinds of currents in the system. 
When the system is isolated from the external circuit, an induced current ac- 
cumulates an extra charge (or spin, orbital) on one side of the system, which 
results in a spatial polarization of charge, spin and orbital. 39 Accordingly, 
this type of current associated with such internal degrees of freedom is often 
called a polarization current. In the following, we describe the physics of po- 
larization current using the language of non-abelian gauge field, in particular, 
that of Tki^^t- One of the advantges of taking such a viewpoint is that the role of 
symmetry becomes transparent, which we summarized as a set of cancellation 
rules in Table 2. 



5.3.1 Charge polarization current 

Let us first consider a charge polarization current. In the case of time- 
dependent perturbation ((HUJl, the EOM analogous to Eqs. (|61l62f) . are found to 
be 



dx^i^^ _ 5eioc(k) 
dt dkfj_ 

dk\l^ 



z 



(^■)t^ (k,i)zO), (81) 



dt 



= 0, (82) 



where j — 1, - ■ ■ ,N . Collecting contributions from all N pseudospin states and 
from all filled k points, one can calculate the charge current carried by the 
ground state as 

/-rpoi\ _ _ f dii \ ^ dx^J"^ 
^ ' " Vbz(2-)^^ dt 
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BZ 



dk 

(2^ 



Tr 



(83) 



where the k- integral was performed over the whole Brillouin zone (BZ). Eq. 
(|83(l is analogous to Eq. (|64|l . which we found for the charge Hall current. We 
can see that the trace of different types of reciprocal field strength, i.e., Tkf,t 
and J-k^kv are related to different types of physical currents, i.e., polarization 
and Hall type currents. 

We have seen in the previous section that no charge Hall current flows when- 
ever the sysytem is invariant under time reversal T. In contrast, we are going 
to see below that the charge polarization current vanishes whenever the system 
is invariant under spatial inversion /, ^® 



Hioc (x, P, t) = Hioc (-X, -p, t) . 



(84) 



In this case, its TV-fold degenerate Bloch functions at k is related to those at 
— k upto a certain SU{N) gauge transformation /i(k,i). 



N 



( - X, a|ui(-k,t)) = ^ (x,a\uj{k,t))hji{k,t). 



(85) 



Since the field strength Tk t is related through Eqs. H7I8|I to those wavefunc- 
tions, Eq. H85|l reduces to the following identity. 



(86) 



Consequently, 



Tr 



-Tr 



•^fc.t(k,i) 



(87) 



Eqs. (|83ll and (|87|l indicate that the charge polarization current always vanishes 
in /-invariant systems. 



5.3.2 Parity polarization current 

We have already encountered a similar situation in the previous section. 
Under the time reversal T, ^k^^k^O^) is transformed to ~- (J-'k^k,A~^)) upto a 
SU (N) gauge degree of freedom, as can seen in Eq. (|66|) . As a result, the charge 
Hall current vanished in T-invariant systems. On the other hand, a spin Hall 
current was robust against T-invariance. The reason was that not only Tk^k^ 
but also the spin operator are odd under the time-reversal, as given respectively 
in Eqs. (jHEl) and IHHll- 

Following the same type of logic, we can expect that an orbital polarization 
current may remain finite irrespective of the /-invariance of -ffioc, as far as the 
associated orbital operator H(x, p) changes its sign under the spatial inversion 

-^^i.e., the underlying crystal structure has centro-symmetric points. 
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I. Accordingly we dub this type of orbital current a parity polarization 
current. 

Expecting that the above analogy is indeed a sensible one, let us further 
analyze the parity polarization current carried by the ground state. Since we 
have defined this orbital current as Eq. H52I53|) . we have to consider an EOM 
for 

= (n(x, p)) - i (nx^ + x^n) . (88) 

We derive their EOM in terms of the decomposition (|33|l . Our local Hamiltonian 

(2) 

is time-dependent, and so is its local Bloch function. Therefore, f2Ln appearing 
in Eq. (|33|l remains finite in general: 

dt ""on+^^on 

N f r 1 

= zV dkdk'alik,t) {Hioc)-iVt,{0]^) (k,k')a„(k',i). (89) 

— / L J mn 

m.n—1 

This equation is analogous to Eq. (|72|l . The covariant derivative Vfc = (x^) in 
Eq. H72|) was replaced in Eq. (|89l) by another covariant derivative w.r.t. time, 
i.e., Vt: 

f*^*) = ^Jr,^n - l( At) . (90) 

V / mn Ot V / mn 

Let us further develop the analogy between the two cases, i.e., we rewrite Eq. 
H89(l in the following way, precisely as we rewrote (|72|l as 1)73(1 . The details of 
the derivation is given in Appendix D, which is in parallel with Appendix C, 
and the result is 

^ = f: /dkp(k)z:,(k,^)(n)_(k)^^.„(k,i) 

m,n— 1 

- iz;„(k,i)((n)(k).Ffe,,t(k) + h.c.) z„(k,i). (91) 

Following the prescription 1)18(1 , we see that the parity polarization current car- 
ried by k Bloch electron occupying the j-th pseudospin state is given by. 



dOfi^^'' (7)t/TT\/'^^ (7) ^^ioc(k) 

= zWt(n)(k)zW 



- (z(^)t(n)(k).Ffe,,t(k)z(^)+c.c.). (92) 

After taking its summation over N pseudo-spin states and over filled k points, 
we finally obtain a parity polarization current carried by the ground state. 



/-[pou _ [ dk ^ dO^ 



BZ (27r)^ f^^ dt 



'See Eq. 
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BZ 



dk 

(2^ 



Tr 



(n)(k)^,,,t(k) 



(93) 



The first term of Eq. (|92|l did not contribute to Eq. H93() due to a cancellation 
associated with the T-invariance, 



BZ 



dk 



Tr 



(n)(k) 



^eioc(k) 



0. 



(96) 



On the contrary, the Berry phase contribution, i.e., Eq. 1931) turns to be quite 
robust against both / and T symmetries. In particular, in the /-invariant case, 
(n)(k) is identical to — (n)(— k) up to the SU{N) gauge transformation h, 

im-k) = -h-\U){k)h (97) 

This can be shown explicitly using Eqs. H5()I57I85I54|I . Eqs. (|86I97|) indicate 



Tr 



(n)(k).Ffe,.t(k) 



Tr 



{n)i-k)Tk,A-^) 



(98) 



Eq. H98(l holds also true in the T invariant case, as is clear from Eqs. 
Eqs. H93I98|I confirm our hypotheses that the parity polarization current is 
indeed robust against / symmetry. 



5.3.3 Quantum spin pump 

Another possible direction to be explored is to study how to induce a spin 
polarization current by breaking both T-invariance and /-invariance. This sce- 
nario can be implemented |15| in a certain kind of quantum spin chains such 
as Cu-bensoate and Yb4As3. The ground state of these quantum magnets is 
known to be quantum critical point (QCP), which is interpreted as a Dirac 
monopole, i.e., a source of the U(l) field strength Tkt- When this quantum 
system is driven around this QCP by applying an electric field E and/or mag- 
netic field B, a spin polarization current can be induced. The electromagnetic 
fields break both the /-invariance and T-invariance. They also induce a spin 
gap and realize a quantum critical point at the origin of E-B plane. When the 
system goes adiabatically around this origin, a quantized number of spins will 
be transported from one edge to the other through the system. This quantized 
value is a physical manifestation of the first Chern number associated with the 
QCP. 

^"T-invariance relates A^-fold degenerate Bloch functions at k with the ones at — k upto a 
SU{N) gauge degree of freedom g as Eq. I65i . This impUes, 

((nX-k))' = g~i(n>(k)3, (94) 

(^fc^t(-k)) = g'^:Fk^t{^)9. (95) 
Eq. I94i justifies 1961 . whereas multiplying Eq. I94i with Eq. I95i one finds immediately 
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5.4 J^fe,iXp associated with the spatial inhomogeneity 

Contrary to J^k^k^ and Tk^t, the reciprocal field strength Tk^x^ does not 
seem to be related directly to a physical observable such as Hall type currents 
and polarization currents. However, when the system contains spatial inhomo- 
geneity such as lattice defects, Tk^x^ appears and plays an important role in 
the dynamics of Bloch electrons around the defects j7j . Another possible appli- 
cation of Tk^x^ is the electron transport properties around a magnetic domain 
wall, where the spatial modulation of ferromagnetic moments induce Tk x^-, and 
naturally influences the EOM for the electron wave packet through this Berry 
curvature. 

6 Conclusions 

We have derived and analyzed the semiclassical EOM for a wave packet of 
Bloch electrons, under perturbations slowly varying in space and in time. Their 
interpretation in terms of non-abelian gauge field in the reciprocal parameter 
space was the central issue of the paper. The same type of EOM has been 
previously derived for the abelian, i.e., U(l) case, by using either (i) Time- 
dependent variational principle ^ or (ii) Path-integral method using Wannier 
basis 1201 ■ We have generalized such EOM to a non-abelian case by using only 
the most fundamental principles of quantum mechanics. 

The advantage of our formalism was that 

1 . it was asymptotically exact in the framework of linear response theory, as 
a result of systematic expansion w.r.t. the perturbation /3 or x — x, 

2. it revelaed that there are different types of gauge field of different physical 
origin, 

3. it was useful for developing symmetry analyses on various types of Berry 
phase transport. 

The first point refers to Eq. (|20|) and all the related analyses developped in 
Sees. 3 and 4. The relevance of our results in relation to other approaches was 
further discussed in Sec. 4.4. As for the second point, two different sources of 
gauge field have been revealed, i.e., (i) Projection onto a subspace spanned by 
N Bloch bands, (ii) Bloch basis moving in the course of time. The former is the 
origin of Tk^k^ which is directly related to spontaneous Hall currents of various 
degrees of freedom. The latter brings about Tk^^x^ and J^fc^t, which plays an 
important role in the spatially and temporally inhomogeneous system. 

Finally, concerning the last point in the above list, we have applied our 
formalism to the analyses on the spin and orbital transport phenomena with 
the help of Boltzmann transport theory. The role of time reversal and space 
inversion symmetries in the appearence of finite Hall/polarization currents has 
been extensively studied. The cancellation rules are summarized in Tables 1 and 
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2. The concept of parity polarization current has also been introduced, which 
may concretize Berry phase engineering in the context of orbital transport. 

We leave for a future study further investigations on their application to the 
domain wall physics and that of quantum pumping. In conclusion, we believe 
that our analyses on non-Abelian gauge field will see in the near future a possible 
application in the context of Berry phase engineering. 

Note added - After completion of this work, we were informed of a related 
effort by D. Culcer, Y. Yao and Q. Niu. gO] 
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Appendix A : Matrix element {H)^^(k) and Afi^^j 

Matrix elements of the linearlized Hamiltonian, i.e., Eqs. (|36I37|I . have been 
extensively used in Sec. 4.3. Let us recall those equations together with the 
matrix elements of noncommutative coordinates, i.e., Eqs. (|23I27I25I28|) . Our 
purpose here is to substitute the expression H37|l into 



N 

E 

771,71 — 1 



dka;,(k) (iJ),(.T^) (k)a„(k), 



(99) 



an expression analogous to the second line of Eqs. (|33|l . and to rewrite 57^^'' in 
terms of the field strength fjFt h 1 

The first term of Eq. H37(l . i.e., eioc(k, x, i) gives a standard velocity term 
when inserted into Eq. H99|) . Since [.F^^fc^] is related to the commutator, 
[Vfc^,Vfe„]^^ or equivalently, [(x^), (a;^)]^^(k), one can easily imagine that 
the last two terms give in Eq. H37|l give when inserted into Eq. ()99|l a contri- 
bution related to \!Fk h 1 • One can indeed verify 



r 1 / 9eioc 



2 V dx^ 
9eioc 



d2 



N 

1=1 



9eio, 



de 



loc 



N 

E 

1=1 







rnl 





In 



dx„ 
i 

+ 2 



In 



dkfj^dxu 



' ku, 



mn dx, 
Vfe. 



(100) 
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The second term of Eq. H37II gives, when inserted into the commutator, 



geioc(k,x, ^) 
dx,, 



dk^dxy 



(101) 



CoUecting the contribution (|101|l and the last two terms of Eq. H1UU|) . i.e., terms 
not related to ^fc„fc^, one defines ^^^x} introduced in Eq. H38|l ; 



dkudxu' 



p(k,i) a;,(i)-«^™(k,i) Vfe„ z„(k,i)L (102) 



This term vanishes after fc-integration with the help of prescription given in Eq. 

m- 

Appendix B : Derivation of the EOM for z(t) 

We demonstrate here the derivation of EOM for z(t), i.e., EOM describing 
the motion of the internal pseudospin degree of freedom. For that purpose we 
once have to go back to Eq. (|30|l . After multiplying it with a weight p(k, t), 
we integrate it over all the k-points, to find. 



/rfk|l|(p(k,t))z(k,t)+p| 
i j (ik/9(k, I - <^off + 

2 1 ''^p^'^^'^^^dk;^ 



■A, 



dXfj 



-Ax,. + A 



(103) 



where eoff — eiocl + Ae as given in Eq. (|^ . In Eq. (|103|l . repeated /^-indices 
should be summed over /i = 1, • • • , £). In order to obtain Eq. (|103|l . we also 
used. 



/■ r aeioc(k) 9 / \ 1 9eioc 1 1 



= 7; d^P 



deioc dz 
dxfj_ dkf^ 



(104) 



We now substitute. 



-(p(k,t))=^— ^p(k,t) 



dku 



^ c)ei oc(k,x, ^) 

dXa 



into Eq. H1Q3|) . then perform a partial integral w.r.t. /c^. The result is, 

deioc dz dz, 



dkp(k, t) 



dx^ dkfj_ dt 



deioc deioc 



= ^ / dkpi -e,f, + x,^- ^A,^ + ^-J^A-x^ + At \z. (105) 



dx^i 
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Finally, in order to rewrite Eq. (|105|l in the form of Eq. (|44f) and complete its 
derivation, we adopt the prescritpion (|18|l . 



Appendix C : EOM for and spin Hall current 

Our purpose here is to rewrite an EOM for Eq. (|72|l into its final form, i.e., 
Eq. H73|l . so that we can express the spin Hall current as the following trace in 
the pseudospin space, Tr[(5Q)(k)J^j,^j_^(k)]i?,y. 

Let us first recall the assumption we made in Eq. H58|l . This assumption 
allows us to factorize (C'^„)^^^^Jk, k') in Eq. (|72|l into a product of {Sa} and 

N 

(S'q.X^ + Xfj^Sa") — ^ ^ (^{S a) ml p,) In + p,) ml {S a) In^ ■ 



1=1 



Correspondingly, the commutator appearing in Eq. (|72|l can be decoupled into 
the following two types of commutators. 



(106) 



The first term together with its hermitian conjugate conceives a commutator 
between (H) and {Sa), which consitutes the EOM for the spin : Firstly, 
we show that this commutator vanishes. Since (Sa) is diagonal in k-space, it 
clearly commutes with (Hq); 



(ifo) , {Sa) (k', k) = 5(k' - k) eioc(k), {Sa) (k) 

J rnn L 

= 0. 



(107) 



Therefore, the commutator between {H) and {Sa) becomes proportional to the 
covariant derivative of {Sa) (k) w.r.t. fc^, i.e.. 



u 

{H),{Sa)\ (k',k) = eVi?J(a;,,),(5„)l (k',k) 

- ran L J rnn 



<5(k'-k)^£;^[Vfc^,(5„) (k) 



(108) 



Because the spin operator itself does not depend on the crystal momentum, the 
partial derivative of (k) w.r.t. fc^ reduces to the commutator between iAu 
and (S„) (k), : 



dk,, 



Sq 



1=1 



Ul){Ul 



Sa 



"n(k)^ + <^Um(k) 
Sa 



Sa 



dk„, I 



Ul){Ui 
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In the second line fc-dependence is not written explicitly. We used Eq. lfSS|l 
between the two lines. Then the covariant derivative ol (Sa) (k) w.r.t. fc^ 
appearing in Eq. H1U8|I also vanishes, 



Vfe,,(5a)(k) 



^( (^o>™n(k)) - ^Av,(5.)(k) 
0. 



(110) 



Consequently, the first term and its hermitian conjugate in Eq. Hl()6|l are indeed 
zero. On the other hand, the second term in Eq. (|106|l contains the commutator 
between (H) and (x^), which describes the EOM for Xi, now. This term gives rise 
to a field strength J-^ through the commutator between covariant derivatives 
w.r.t. different components of the crystal momentum, i.e., [Vfc^,Vfc^]. Namely, 



{H) , {x.) 



D 

(Ho) +eJ2Et.{ 

D 



(k',k) 



(5(k'-k) [eioc(k)l + e^£;^ zVfc^,iV;. 
-^50^ - k) j^l - ef] £;,(.F..,...„(k)) |. (Ill) 



Finally, substituting Eq. ((TU5|l together with Eqs. (|108I110|I and ((TTT|l into Eq. 
(O, one finds Eq. GSJ. 

Appendix D : EOM for and parity polarization current 

In parallel with Appendix C, we rewrite below Eq. H89|) into its final form, 
i.e., Eq. H91|l . Let first recall the assumption (|59|l . which says that the parity 
operator 11 has no matrix element outside the iV-fold degenerate band. This 
implies that (Oj?) can be factorized into a product of two N hy N matrices, 
or (na:^) = (n)(; 

can be decomposed into two types of commutators as 



Thanks to this factorization, the commutator in Eq. H89|l 



+ i((n>[(ffioc>-iVt,(a;^)] +h.c.) (112) 



On the r.h.s., the first term is a commutator between (i/ioc) — i^t and (11), 
which constitutes the EOM of the parity under time-dependent perturbations: 
dH/dt. Since the parity operator itself is independent of time, we can prove 
that this commutator indeed vanishes, in the same way as we did in Eq. Ij 109(1 . 



[(i/ioc) -*Vt,(n)] =0. 



(113) 



On the other hand, the second line of Eq. H112|l is a commutator between 
(Ho) — iVf and (Xfj), which describes the EOM for Xfj,. This commutator gives 
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rise to T^^t through a commutation relation between two covariant derivatives, 
one w.r.t. time and the other w.r.t. the momentum, [iV/t^,zVt] =i^k^t- Thus 
the second hne of Eq. 1112|) may be rewritten as 



(k', k) = iSik' - k) eioc(k) - iVt, i^k. 



= -i5(k' - k) I 



9eioc(k)^ 



(114) 



Substituting Eq. ifTHIl together with Eqs. (|113I114() into Eq. one finds 

Eq. jnH)- 
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invariance 


T 




charge : 1 


+ 


+ 


spin : (5) (k) 




+ 


parity : (n)(k) 


+ 




Hall-type current: Tk^k^O^) 




+ 


polarization current: J-k^tO^) 


+ 





Table 1: Transformation properties under time reversal T and spatial inversion 
/ - A negative (positive) sign indicates whether a matrix element in question 
at k reverses its sign (or not) in comparison with that of — k when the system 
is invariant under a certain symmetry operation, such as T or /. For example, 
the negative sign for (5*) (k) in the T- invariant case means that the spin matrix 
(5*0) (k) is identical to — ( (Sa) (— k))* upto a certain SU{N) gauge transforma- 
tion as given in Ea. l|t)8|) . 



type of current 


Hall-type 


polarization 


invariance 


T 




T 




charge 




+ 


+ 




spin 


+ 


+ 






parity 






+ 


+ 



Table 2: Cancellation rules for charge/spin/parity Hall-type/polarization cur- 
rents - A negative (positive) sign indicates that contributions to the total 
charge/spin/parity (vertical axis) Hall- type/polarization (horizontal axis) cur- 
rent from k and — k electrons (do not) cancel each other when the system is 
invariant under either T or /. This table can be deduced from Table 1. For 
example, in the T-invariant case, a negative sign for the spin in Table 1 gives, 
together with another negative sign for the Hall- type current, a positive sign for 
the spin Hall current in Table 2, correspoinding, respectively, to Eqs. H68I66|I 
and 
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